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Abstract. We prove a two-term asymptotic expansion of eigenvalue sums of the Lapla- 
cian on a bounded domain with Neumann, or more generally, Robin boundary conditions. 
We formulate and prove the asymptotics in terms of semi-classical analysis. In this refor- 
mulation it is natural to allow the function describing the boundary conditions to depend 
on the semi-classical parameter and we identify and analyze three different regimes for this 
dependence. 



1. Introduction and main result 

1.1. Introduction. The Laplace operator on a bounded domain $7 C M. d , d > 2, initially 
defined as a symmetric operator in L 2 (f2) with domain C^(Q), admits various self-adjoint 
extensions that correspond to different boundary conditions. Our goal in this paper is to 
study how different boundary conditions influence the asymptotic behavior of the eigenval- 
ues. 

We consider self-adjoint extensions that are generated by a quadratic form 

[ \Vv\ 2 dx + [ c(x)\v(x)\ 2 da(x), ceij'ffi). (1.1) 
Jn Jan 

Here the form domain H 1 ^) is the Sobolev space of order 1, da denotes the d— 1-dimensional 
surface measure on the boundary dQ, and c is a bounded, real valued function on dQ. This 
quadratic form induces a unique self-adjoint operator — A c in L 2 (fl) and functions from the 
domain of — A c satisfy, in an appropriate sense, Robin boundary conditions 

dv 

—— (x) =c(x)v(x), xedtt, (1.2) 
on x 

where denotes the inner normal derivative. We remark that c = corresponds to the 
important case of Neumann boundary conditions. The Dirichlet Laplacian, generated by the 
quadratic form \Vv\ 2 dx with form domain Hq(Q), can be recovered formally by taking 
the limit c — > oo. 

If the boundary of f2 is sufficiently regulary (e.g., Lipschitz continuous), the spectrum 
of — A c is purely discrete: It consists of a sequence of eigenvalues Ai < A2 < A3 < . . . 
that accumulate at infinity only. Here we study how the asymptotic distribution of the 
eigenvalues depends on the boundary condition induced by the function c. 
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It is a classical result that the eigenvalues satisfy 

^ = J^ W - d n^ + o(n^) as n ^ oo , (1.3) 

where |0| is the volume of and denotes the volume of the unit ball in M. d . In the case 
of Dirichlet boundary conditions these asymptotics go back to |Weyl2| . They have been 
generalized in various ways, in particular, to the case of Robin boundary conditions (II. 2j) ; 
see, for instance, the lecture notes [BS80J. 

It has been conjectured by Weyl that (|1.3|) is the beginning of an asymptotic expansion in 
n and that the second term should depend on the surface area of Q. Initially, a weaker form 
of this conjecture has been verified, not for individual eigenvalues, but for smooth functions 
of the eigenvalues; see, e.g., [Ple54l lMS67]. For instance, |BG90] computed in the case of 
boundary conditions (|1.2|> 



V e _tA * = (47rt)~ d/2 ( \tt\ + ^\dn\ i 1 / 2 + \ I (H(x) - 6c(x)) da{x) t + 0(t 3/2 )) 
~{ V 2 3 Jan J 

as t->0. (1.4) 

Here H(x) is the mean curvature (the trace of the second fundamental form) at x £ d£l. 
We see that the second term indeed depends on the surface area \dd\ and is independent of 
c. The boundary condition enters only in the third order term. (For Dirichlet conditions, 
however, the sign of the second term flips.) In contrast to (11. 3p . the expansion (|1.4p requires 
the boundary to be smooth. 

A two-term asymptotic formula for individual eigenvalues was eventually shown in a cel- 
ebrated work of V. Ivrii; see |Ivr80allFr^^lSV971IIvr98] . He showed that, under a certain 
condition on the global geometry of $7 (and some smoothness conditions), one has for bound- 
ary conditions (II. 2ft 



A - 4?r n 2/d 271 n i/d | o(n i/d, ag _ .„ fl5l 

Again, for any bounded function c the result is the same as for Neumann conditions. We 
emphasize that (11.5P implies the two-term analogue of (jl.4p . but not vice versa. 

In this paper we shall study an eigenvalue quantity which is intermediate between (jl.4|) 
and (|1.5p . namely, partial sums Y^=i^j as n — > oo or, equivalently, X^=iO\? — A 4 )- as 
fi — > oo. These partial sums describe the energy of non-interacting fermionic particles in 

at fixed particle number n or at fixed chemical potential fj,, respectively. They play an 
important role in physical applications. 

Since the function A i-> (A — fx)- is not smooth, we cannot expect that a three-term 
asymptotic expansion exists for these eigenvalue sums. Hence, to see the effect of boundary 
conditions already in the second term of the asymptotic expansion we have to choose energy- 
dependent boundary conditions. Let us state this problem in a semi-classical set-up. For 
a small parameter h > we define self-adjoint operators H(b) = —h 2 A b / h — 1 in L 2 (Q) 
generated by the quadratic form 

q b [v] = h 2 [ \Vv\ 2 dx + h [ b(x)\v(x)\ 2 da(x) - [ \v(x)\ 2 dx (1.6) 
Jn Jen Jn 
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with form domain H l (£l). Here b is a bounded function on dfl that may also depend on h. 
The quadratic form q b induces, in an appropriate sense, /i-dependent boundary conditions 

h-p-(x) = b(x)v(x), x£dn. (1.7) 

In this introduction, we denote by E n (b, h) the eigenvalues of the operator —h 2 A b / h ; conse- 
quently, the eigenvalues of H (b) are given by E n (b, h) — 1. As we explained, our main goal 
will be to study the sum of the negative eigenvalues of H(b), 



TrH(b)_ = J2( E n(b,h)-l) 



neN 

in the semiclassical limit h 1 0. We prove two-term asymptotics and show how the second 
term depends on the function b. Our analysis will show that the asymptotics has different 
forms in three different regimes depending on the size of b as h \, 0. The three different 
regimes are where b — > as h 1 0, b of order one as h I and \b\ — > oo as h I 0. 

As an example of the first regime, let us consider the case b = he with a bounded function 
c independent of h. This corresponds to the classical situation discussed above, where the 
boundary condition (II. 7ft is independent of h and therefore the eigenvalues E n (b, h) = h 2 \ n 
depend trivially on h. Then fll .5f) implies 



iE^=(^^»--^(5f^^"-M»''«) - ^00,(1.8, 



and this is equivalent, by a simple majorization argument, to 

TrH(b)^=L^ ) \n\h- d + ^L i jl 1 \dn\h- d+1 + o(h- d+1 ) as | (1.9) 

with = ^2 (2-7r) _c! ci;rf. Of course, we find again that the first two terms of the asymptotics 
are independent of the boundary condition. As we shall see, this is characteristic for the 
whole regime where b — > as h \. 0. We emphasize that as a byproduct of our analysis 
we establish (II. 9p independently, without using (jl.Sp : see Theorem 11.21 This includes, as a 
special case, the Neumann Laplacian. 

Among the three regimes mentioned above, the technically most interesting one is when b 
is independent of h. In this case the second term of the semi-classical limit of TrH(b)- does 
depend on the local behavior of b(x); see Theorem 11.11 b elow . 

Finally, in Theorem II .3} we consider functions b such that \b\ diverges as h \. 0. In this 
case, the form of the asymptotics depends on whether b is negative somewhere or whether 
b is non-negative. In the first case, the asymptotics are determined by the negative part of 
b alone. Moreover, if b diverges fast enough, then the boundary term becomes the leading 
term and diverges faster than the Weyl term. On the other hand, when b is non-negative 
the order of the second term is preserved but the coefficient may change. 

We obtain these results by further extending the approach developed in [FGll|,lFG12] . 
where we treated the Dirichlet Laplacian and the fractional Laplacian on a domain. One 
virtue of this approach is that it requires only rather weak regularity assumptions on 30 
and b. Essentially, a C 1 assumption on dQ and on b suffices for a two-term asymptotics. 

We now turn to a more precise description of our assumptions and results. 
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1.2. Main Results. Let C Mr, d > 2, be a bounded domain such that the boundary 
satisfies a uniform C 1 condition. That is, the local charts of dQ are differentiable and their 
derivatives are uniformly continuous and share a common modulus of continuity; see (14. lh 
for a precise definition. Moreover, we assume that the boundary coefficient b is a continuous, 
real-valued function on <9J7 and we denote a modulus of continuity by (3, i.e., 

\b(x)-b(y)\<P(\x-y\) (1.10) 

for all x, y £ <9$7. We assume that /3 is non-decreasing. 

We remark that the boundary conditions fjl -Tj) for functions in the operator domain of 
H(b) need not hold in the classical sense under these weak assumptions on the boundary. 
For dQ € C , however, this operator can still be defined by means of the quadratic form 
and the characterization of the operator domain in terms of the form domain gives a weak 
sense in which (ll.7p are valid. This suffices for our proof. 

For a constant b E M. we set 

(C d (-f + Jj(l - P 2 )( d+1 )/ 2 ^p) for b > , 

4 2) ( & ) = |c d f for 6 = 0, (1.11) 

k Q (-f + ^(l - p 2 )^ 2 ^^ + vr(6 2 + l)( rf +D/ 2 ) for b < , 

where = 4\§> d ~ 2 \(2ir)~ d (d 2 — 1) —1 . This expression comes from the explicit diagonalization 
of a one-dimensional model operator; see Section [3l Although it is not obvious from the def- 
inition, the function L d '(b) is continuously differentiable and non-increasing; see Lemma [3.5l 
and the remark after Proposition 13.11 In particular, for b > 0, we have 

" = i™ L ?\b) < L d 2 \b) < limLf (6) = 4 2 )(0) = (1.12) 

with defined after (fL9|) . 

To control error terms we have to introduce a non-decreasing function 5 : [0, ||6||oo] - > 
[0, oo) such that 

5{X) >\{xedn : 0< \b(x)\ < X}\ (1.13) 

for all < A < ||6||oo- 

Our first main result is the following. 

Theorem 1.1. Let dfl G C 1 and assume that b satisfies (jl.lOj) and (|1.13p with (3(1) = o(l) 
and 5(1) = o(l) as 1 1 0. We write 

Tr(#(6))_ = 4 1} 1^1 ^ d + I Lf (b(x))da(x) h~ d+1 + R h . 

J dQ. 

Then, for an h-independent domain f2, a given h-independent upper bound on ||6||oo an d 
given h-independent /3 and 8, the asymptotics 

R h = o(h- d+1 ) 

holds uniformly in b satisfying these conditions. 

In other words, in this theorem we claim that Rh = o(h~ d+1 ) if b is independent of h. 
Moreover, we claim that these asymptotics are valid even if b depends on h, as long as it 
can be controlled in some uniform way. More precisely, we prove that given f3 and 5 (both 
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non-decreasing and vanishing at zero) and constants C > and e > 0, there is an h e > 
such that \R h \ < eh- d+1 for all < h < h £ and all b satisfying \\b\\oo < C, (fl"T0|) and (fL13]) . 
Our proof would also allow us to consider /i-dependent domains but we do not track the 
dependence of the constants in terms of 0, for the sake of simplicity. 

Our next result concerns the case where ||6||oo as /i | 0. We will see that the 
asymptotics are the same as in Theorem 11.11 with 6 = 0. We cannot apply Theorem 11.11 
however, since for 6^0we cannot choose 5 independent of h such that (II. 13ft is satisfied and 
5(X) = o(l) as A J, 0. Moreover, we can dispense with the assumption that b is continuous. 

Theorem 1.2. Let <9$7 £ C 1 and assume that b = 6(h)bo with 6(h) = o(l) as h \, and with 
a bounded function bo . We write 

Tr(tf(6))_ = Lf |fi| hT d + ^4-1 W h ~ d+1 + R h ■ 

Then, for an h-independent domain f2 and a given h-independent upper bound on ||&o||oo> 
the asymptotics 

R h = o(h- d+1 ) 
holds uniformly in b satisfying these conditions. 

We refer to (|2.17p for an explicit bound on R^. 

Our third result concerns the case where b = O(h)bo with Q(h) —> oo. 

Theorem 1.3. Let dQ E C 1 . Assume b = @(h)bo with 0~ 1 (/i) = o(l) as h J, and with b^ 
satisfying (jl.lOp with (3(1) = o(l) as I \, 0. We write 

Tr(iT(6))_ = 4 1} 1^1 h ' d + TrCd / + R h . 

Jdtt 

Then, for an h-independent domain Q, a given h-independent upper bound on ||6q||oo and a 
given h-independent f3, the asymptotics 

R h = (Q(h) d+l h- d+l ) 

holds uniform in b satisfying these conditions. 

If, in addition, b(x) > for all x G dfl, 5(\) = o(l) as X 10, and j3(Mh)Q(h) = o(l) as 
h\.Q for every fixed M > 0, then 

Tr(F(6))_ = 4 1} \V\ h~ d - -4!^i h- d+1 + ^4-i l^ol h- d+1 + o(h- d+1 ) , 

where dVt + = {x e dfl : b(x) > 0} and dflo = {x G dfl : b(x) = 0}. 

We emphasize that, if the negative part of b does not vanish and 0(h) = h 1 with 7 = 
l/(d + 1), then the order of the boundary term is the same as the order of the Weyl term. 
For 7 > l/(d + 1) the boundary term becomes the leading term. 

Since j3(l) vanishes at most linearly in I for non-constant 6, the condition f3(Mh)0(h) = 
o(l) as /i ^ in the second part of the theorem implies 0(h) = o(/i _1 ). Our techniques do 
not allow us to consider faster growing 6's and we do not know whether one still can expect 
the result in that case. 
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2. Strategy of the proof 

In this section we outline the main steps of our proof. In particular, we explain how the 
main results follow from local estimates. 

First, we localize the operator H(b) into balls, whose size varies depending on the distance 
to the complement of $7 [H6r85, SS03]. Then we analyze the local asymptotics separately in 
the bulk and close to the boundary. 

To localize, let d{u) = inf{|x — u\ : x ^ fi} denote the distance of u £ M. d to the 
complement of f2. We set 

i( u ) = l -(i + (d{uf + iiy 1/2 y\ (2.i) 

where < Iq < 1 is a parameter depending only on h. Eventually, we will choose Iq = o(l) 
as h \ 0. In Section [5] we introduce real- valued functions 4> u £ Co°(M d ) with support in 
B u = {x £M. d : \x — u\ < l(u)}. For all u 6 M. d these functions satisfy 

Hu\L<C, IIV^I^ KCliu)- 1 (2.2) 

and, for all xeK d , 

/ (p 2 u (x)l(uy d du = 1. (2.3) 

Here and in the following the letter C denotes various positive constants that are independent 
of n, Iq and h, but may vary from line to line. To estimate error terms in the following results 
we put 

b m = inf b(x) . 
x&an 

Proposition 2.1. There is a constant Cq > such that for < Iq < C^ 1 and < h < Iq/4 
the estimates 

-C(l + (b^hlQ 1 ) l^h-^ 2 < f Tr {<j )u H{b)<t> u )_ l{u)' d du - Tr(H(b))_ < 

hold. 

This proposition will be proved in Section [5j 

In view of this result one can analyze the asymptotic behavior of Tr(0^^(6)(/> u )_ sep- 
arately on different parts of S7. First, we consider the bulk, where the influence of the 
boundary is not felt. 

Proposition 2.2. Let 4> £ Cq(Q) be supported in a ball of radius I > and let 

IIV0IU < c^r 1 . (2.4) 

Then for all h > the estimates 

< L { j } [ (j) 2 {x)dx h~ d - Tr {4>H{b)<j>)_ < Cl d ~ 2 h- d+2 (2.5) 
Jn 

hold, with a constant C > depending only on CV 

For <j> G Cq(Q) we have (j)H{b)4> = <p(—h 2 A — l)c/>, where —A is defined on the whole 
space L 2 (M. d ) with form domain ff 1 (IR d ). Hence, this result is independent of the boundary 
coefficient b and the proof of Proposition 12.21 is the same as in [FGllj . 
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Close to the boundary of Q, more precisely, if the support of (j) intersects the boundary, a 
term of order h~ d+l appears that depends on 6. In this situation let B be a ball containing 
the support of and put 

b-= inf b(x), V= inf \b(x)\, b s = sup |6(x)|. (2.6) 
xednnB xednns xednnB 

To state the remainder estimate we denote by w a modulus of continuity of the boundary of 
f2; see (|4.ip for a precise definition. 

Proposition 2.3. Let cj) E CQ(JH d ) be supported in a ball of radius I > and let inequalities 
(jl.lOP and (|2.4p be satisfied. Then there is a constant Cq > such that for < I < C^ 1 
and < h < I we have 

,,, M , /J( , ; L ril /,.-'( ,-w7,.;-,-''j_ / f (2) iJ,i,')i^ir'i,^i,w>-'-'- 1 

with 



Tr (cf>H(b)(f>)_ =Ly / <j)\x)dxhr d + / L]; , (b(x)) ( l) 2 (x)da(x)h- d+1 +R bd (h,l,b-,b t ) (2.7) 
n J9o 



\d+l 



\R bd (M,b-,v)\ < C ^ d +o,(z) (i+^ ( 

+j(l + (6-)l)^(0 
For b s < /i/Z we a/so /love 

Tr(0if(6)0)_ = 4 1} / cp 2 {x)dxh~ d + 7^-1 / </> 2 (x)aV(x) + i? (/t, /, b s ) (2.8) 

|i? (/i,^,o s )| < C7 d /T d (r 2 /i 2 + 6 S (1 + I In 6 s I ) • 
Here the constants C > depend only on VI and CU. 

The first statement in Proposition 12.31 is the crucial result of this section. It yields a 

(2) 

precise estimate with the boundary term including the correct constant L d (b). However, 
we obtain an error term that diverges as b l — > 0. To overcome this effect we also need the 
second statement for b very close to zero. The next lemma is a simplified version of (|2.8p . 
where we estimate the boundary term by Cl d ~ 1 h~ d ~ >r ^ . 

Lemma 2.4. Under the conditions of Proposition DOl there is a constant Cq > such that 
for < I < 1 and < h < I we have 

Tr(0tf(6)0)_ = 4 1} f 4> 2 (x)dxh~ d + R' {h,l,b~) (2.9) 
Jn 

with 

\R' Q (h,l,b-)\ < Cl d h- d (r x h + u{f) + r 1 h{b-)i +1 (minllh- 1 ^-)-,!} +w(l))) • 

Both Proposition 12.31 and Lemma 12.41 will be proved in Section RJJ 

Based on the preceding results we can now give the proofs of our main results. 

Proof of Theorem \l.l[ We fix two parameters < A < 1 and < (jl < 1/4 and set Iq = hfi^ 1 . 
Let us recall the definition of l(u) from (|2.ip and of B u = {x S M. d : \x — u\ < l(u)}. We set 

u = {u g R d ■. dnnB u ^q>}. 
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First, we need to estimate l(u) uniformly. Note that by definition 
l( u ) — — min 1) and l(u) > — >h 



(2.10) 



for all u G M. d . Moreover, for u G U, we have d(u) < l{u) and 

l(u) < l /V3 = h/(VSfJ.) . (2.11) 

For < h < ^Cq 1 it follows that £0 < C^ 1 an d Z (u) < C^ 1 for all u £ U. Moreover, 
h = uIq < ^o/4 < /(u). Therefore the assumptions of Proposition \2.1\ Proposition 12.21 and 
Proposition 12.31 are satisfied. 

Depending on A we decompose U into the regions 



U = {ueU -. 3xednnB u 
u* = {ueu :Vi£ onnB u 
U > = {u£U : 3 x e an n B u 

We remark that U = Uq U U* U [/> and that the three sets are mutually disjoint. Indeed, if 
x G 30 H 7? u with it G Uq, then by the continuity of b, see (jl.lOp . 

ft 



b(x) = 0} , 

< < A} , 

|6(x)| > A} . 



\b(x)\<P(l(v))<P 
and similarly, if x G <9f2 fl B u with it G C7> , 

|6(x)| > \-p 



\Z3fxJ 



(2.12) 



Thus, by our assumption on /3, we have for all sufficiently small h > (depending on \i and 
A) that /3 (j^jij < X - (3 (7^7)- Thus [/ n[/> = 0, as claimed. We can also make sure 



for all x <E dfl D B u with u G ?7o 



that for all sufficiently small ft 

\b(x)\ < v^/x < ft/Z(it 
and 

\b{x)\ > A/2 for all x £ dtt n B u with it G U > . 
To estimate error terms we put, similarly as in (12. 6|) . 

ft- inf ft(.r). ft', ■ inf |ft(.r)|. ft 



(2.13) 



inf 6(x) , &?, 
x£dnr\B u 



inf 



sup \b(x) 

xednnBu 



First, we apply Proposition 12.11 Then, in order to estimate Ti(<p u H (b)cp u ) - , we use ([2.5 
for u G n \ U, (H2D for u G U>, $ZM) for u G Uq, and ([221) for u G 17*. We obtain 



-i? < L 
with 

R 



(i) 



m d Jn 



2, dxdu 



U> 



+ 



+ 



R bd (h,l(u),b u ,bi) 
Lf{b{x))- 



u J on 



du 
T{u) c 



+ 



U 



\R (h,l(u),b s u )\ 



du 

T{uy 



9 . , da(x) du 



+ 



u* J an 



liuYh 11 - 1 Ju 

lu-d+2 



\R' (h,l(u),b u )\ 



1 + (6 m ) d+1 W - 
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I? 



f \R bd (h,l(u),b-X)\i^+ [ \Ro(h,l(u)X)\ 
Ju > 1 H«r JUo 



du 



+ 



+ 



u Jon 



l(u) d 

Lfm)--Al 



i( u y 



r, , , da(x) du 



l{u) d h d ~ l Ju* 



\R' (h,l(u),b v 



du 



u* Jdn 



Lf\b(x)) <fi(x) 



da(x) du 
l(u) d h d 



-i 



+ C 



n\u 



l(u) 2 duh' 



-d+2 



In the main term we change the order of integration and use the partition of unity property 
(ET31) to obtain 



L 



(i) 



<Pl{x)dx^- d h- d = L^Mh- d 



and 



Thus, we get 



u Jdn 



L { f' \b(x))da(x)K 



-d+l 



-r~ < L d 1) \n\h~ d + 



L d 2) \b{x))da{x)h 



-d+l 



Tr(H(b))_ < i? H 



and to complete the proof it remains to bound the remainder terms R^. 

We now argue that the last term in the definition of R + is controlled by the last term in 
the definition of R~, that is, by 



Cl L h- a+ ' z (1 + (6 m ) d+1 W 1 ) < Ch- d+1 n ( 



1 + 116 



I d+l 



To prove this, we note that for u € O \ U we have d(u) > l(u) > Zo/4 and 

t- 2 \dn t \dt) . 



[2.U) 



[ l{ U y 2 du < C (1+ [ d{u)- 2 du) <C [1+ f 

Jn\u \ J{d{u)>i Q /4} J \ Ji c 



h/4 ) 

Here \d£lt\ denotes the surface area of the boundary of Of = {x S : d(x) > t}. Using the 
fact that | d£lt\ is uniformly bounded and that |(?Qj| = for large t, we get 



l(u)' z du < CIq 1 < CfJLh 



n\u 



-i 



(2.15) 



This proves that the last term in i? + is bounded by (|2.14p . 

To proceed, we note that inequalities (|2.1ip and (|2.10p show that l(u) for u G U is 
comparable with Iq. Since B u n dfl ^ we find d(u) < l(u) < CIq and, for any positive and 
non-decreasing function r, 



r{l{u))du < Cr(Cl ) 



du < Cr(Cl )l 



(2.16) 



'{d(u)<l } 

Thus, if we insert the identity l = ha" 1 and the estimates ([27T2|) . (pT3|) . (f2T6|) and (pTT5j) 
into the expressions for R~ and we find that both are bounded by a constant times 



R = h- d+1 (1 + 



+ h 



-d+l 



I d+l 
I oo 



A* 

/i + — + w 
A 



A* y A* VP 



£H~ ( 1 + w 



M / M/ A* \V3/V 



1 + 



ln/3 
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Here we used the facts that \U \ < \U\ < Cl and \L$\b(x)) - \L^ X \ - Cfi(h/y/3/i) for 
x £ B u t~)dQ with u G U . 

To estimate \U*\ we apply Lemma [A. H given in the appendix, to the set N = {x £ d£l : 
< \b(x)\ < A}. By the defining property ()1.13j) of 5 we obtain 

lim sup £ | U* | = lim sup - \U* \ < C5{\) . 
Hence, by our assumptions on to and /3, it follows that 

hmsup (h^R) < (i + wbw^ 1 ) r^+a+c^A). 



A 



By our assumption on 5, the right hand side can be made arbitrarily small by choosing first 
A small and then /j, small. This completes the proof of Theorem 11.11 □ 

Proof of Theorem \1.2\ This proof is similar to the proof of Theorem 11.11 above. Again we 
choose 

U = {u€R d : dttr\B u ^$}. 

and we assume that 1$ = h/j.^ 1 with < \i < 1/4. Then h < l(u) for all u £ U. 

Let us choose h small enoug h such that \b(x)\ = \b (x)\0(h) < V3^ < h/l(u) for all x G d£l 
and u G U. Then we can apply ()2.8[) to estimate Tr (<p u H (b)<p u ) _ for u G U. This yields 



< / \R (h, J(„), 6£)| ^ + CV/^ 2 . 



1/ 



Tr(il(6))_ - 4 1} |0|^- d - -L^il^l^ 41 
Similarly as above we bound 

/ \R (h,l(u)X)\ 1 f^<Ch- d+1 L + U (^)\ + ||&||oo(l+|MI&lloo|)- > ) • 

We multiply this by h d+1 and let h 4- recalling that ||6||oo = #(^)||Mloo = o(l). Since can 
be chosen arbitrarily small, we obtain the claimed asymptotics. □ 

In this case the proof shows that the remainder from Theorem 11.21 can be estimated 
as follows. For all < n < 1/4 we have 

\Rh\ <Ch~ d+1 L + U (^\l + 9(h)\\b \\oa(l + |ln(0(/ l )||6 O ||oo)|)^) • (2.17) 

Proof of Theorem M.SX First, we assume that the negative part of b does not vanish. Then in 
the same way as in the proof of Theorem 11.11 we fix parameters < A < 1 and < \x < 1/4 
and set Iq = /i/x -1 and 

U = {u G R d : dttn B u ^ 0} . 

Here we choose 

u* = { u e u : ix ednnB u : \b(x)\ < A}. 

Then, similar as in the proof of Theorem II .1} by applying (|2.5[) for u G Q\U, (|2.7|) for 
u G U\U*, and (f2T9|) for u G U*, we obtain 



Tr(#(6))_ - L^p\h- d - f Lf \b{x))dxh- d+l 

Jan 



< CR 
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with 



R =h~ d+1 ( 1 + 



I d+1 
loo 



fj,+ — +u 



fCh\ 1 



- 



+ 



+ K 



-d+l 



i + e(h)p 



(Ch 



d+l 



\ n J fi 1 + 
Ch\ 



@(h) Q fCh 



l + U) 



V J 



We emphasize that in order to arrive at this bound we used the estimates \U*\ < \U\ < CIq 
and 



\Lf{b{x))\<C[l + Q{h)(3 



h 



y/3fi 



d+l 



for x G dO, n B u with u E U* . (Note also that the role of (3 in Proposition 12.31 is now played 
by Q(h)f3.) 

To simplify the main term we note that L^\b) = C d nb d+1 + 0(®(h) d ~ 1 ) as h i 0. Hence, 



Lf{b{x))dx = C d 7T I b(x) d _ + 'da(x)h- d+i + 0(&(h) d - L h- d+ ') 



an 

It remains to note that 



\d+l 



-d+l 



\d-l h -d+l-, 



hlO 



mi 



limsup (h d - 1 e(hy d - i R) <C(fi + ^ 



A 



can be made arbitrarily small. (Since we only assume an /i-independent upper bound on 
ll&olloo, one needs to distinguish here the cases whether liminf _1 (1 + ||6||oo) is positive or 
zero.) 

We now turn to the proof of the second part of the theorem. If the boundary coefficient 
b is non-negative we argue in the same way as in the proof of Theorem 11.11 We obtain 

Tr(H(b))--L l p\n\h- i - [ L d 2 \b(x))da(x)h- d+1 

Jan 

*»*(>*Mt)^(t) + KMt)*) 

In ( 0(/i)/3 



+ ~e(h)p 



1 + 



dty c 
4 2) (b(x))da(x 



In this case the continuity of L K d ^(6), see (|1.12p , implies 



i 

'an 4 4 

by dominated convergence as h 4- 0. Again applying Lemma I A. II in the same way as in the 
proof of Theorem 11. 11 we see that all terms equal o{h~ d+1 ) as h \. 0. □ 



To summarize this section, we have reduced the proof of our main results to the proof of 
Proposition 12.11 Proposition 12.31 and Lemma \2A[ 



3. Local asymptotics in the half-space 

From a technical point of view, this section is the heart of our proof. We analyze in 
great detail a model operator which is explicitly diagonalizable. More precisely, we prove 
local estimates corresponding to Proposition 12.31 in the case where f2 is the half-space Mi = 



12 R. L. FRANK AND L. GEISINGER 

{(x',Xd) G R^ 1 x E + } and the boundary coefficient b does not depend on x. Let H + (b) = 
— h 2 A — 1 be the self-adjoint operator in L 2 (W^) generated by the quadratic form 

q+[ v ] = h 2 I \Vv(x)\ 2 dx + hb I \v{x',0)\ 2 dx' - I \v(x)\ 2 dx 

with form domain i? 1 (IR f |_) and with a real constant b independent of x. 

3.1. Statement of the results. Our goal in this section is to prove the following 

Proposition 3.1. Assume that b G R is constant. Let 4> G Q)(R rf ) be supported in a ball of 
radius I > and let (|2.4p be satisfied. Then for h > 

Ti(<f)H + (b) 0)_ = L<p [ cp 2 {x)dxh- d + 4 2) (6) / 2 (^O)cte'/r d+1 + 14 s (M,6) 

with 

\R hs (h,l,b)\ <Cl d - 2 h- d+2 
For \b\ < h/l < 1 we also have 

Tr (<j>H+(b)<j))_ =4 1} / 4> 2 (x)dxhr d + j4-i / ^ 2 (ar',0)d!c'/i- d+1 + J2{ la (/i J i,6). 



1 + 6' 



|< s (/i,/,6)| < C/ d - 2 /i- d+2 (1 + l 2 h~ 2 \b\{l + |ln|b||)) . 
i/ere the constants C > depend only on d and Ga,. 

(2) 

Remark. The proposition shows, in particular, that L d (6) is non-increasing. Indeed, for 
given boundary coefficients b < b' the variational principle implies Tr(iT(6))_ > Tr(H(b'))- 
for all h>0, and Proposition O thus yields L { ^\b) > L { p{b'). 

The first part of Proposition 13.11 is the key semi-classical estimate that we will later 
generalize to curved boundaries and variable 6's. The problem with this bound, however, is 
the in the error term which blows up for small values of b. For that reason we need to 
include the second part, which deals with small values of b. (In passing, we note that since 
L^\b) is continuously differentiable with (0) = as we will see in Lemma 1331 the 

constant i-^-l m ^ ne secon d P ar t °f Proposition 13.11 can be replaced by L^\b) without 
changing the form of the error term.) 

To deal with the transition region between \b\ > 1 (where the first part of Proposition 13.11 
applies) and |b| < h/l (where the second part applies) we need the following rough estimate. 

Lemma 3.2. Assume that b G R is constant. Let 4> G Cq(R ) be supported in a ball of radius 
I > and let (|2.4p be satisfied. Then for all < h < I we have 

Tr(<f)H + (b)(f))_ = 4 1} / <p 2 (x)dxh- d + R'l s (h,l,b) 

Jm.'l 

with 

\R'l s (h,l,b)\ < Cl^hT** 1 (l + 6 d+1 min{6_//i- 1 ,l}) . 
Here C > depends only on d and Ca- 
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From this lemma we immediately deduce a simple bound that will be useful in the following 
sections. 

Corollary 3.3. Assume that b G R is constant. Let (f> G Cq(W 1 ') be supported in a ball of 
radius I > and let (|2.4p be satisfied. Then for all < h < I the bound 

Tr (6)0) _ <Cl d h- d (l + b^hV 1 

holds with a constant C depending only on d and C^. 

The next remark will be used at several places without explicit mentioning in the proofs 
of Proposition 13.11 and Lemma 13.21 

Remark. When bounding error terms in the following proofs we will sometimes encounter 
the term ||</>||oo, which is not mentioned in Proposition 13.11 and elsewhere. The reason is that 
it can be controlled in terms of Cm. Indeed, for x in the support of eft we can choose y at the 
boundary of the support with \x — y\ < I and use (I2.4p to estimate 

|0(x)| = |#r) - 4>(y)\ < \\VHoo\x -y\<C <t> . 
Hence, ||0||oo < C^, as claimed. 

3.2. Analysis of a model operator on the half-line. The bounds in Proposition 13.11 
and Lemma 13.21 are based on the following results about the one dimensional operator — ^ 
on the half-line with boundary condition 

d t v(0) = bv(0) , be R. (3.1) 

For t > and b G M we define 

^(t) = 7r L=cos( i ) + 7 =L=sin(t) 

and, for b < 0, 

$ 6 (i) = v /z 26e 6t . 

In order to treat positive and negative b without distinction we set = for b > 0. Then 
we have 

-d 2 Mt) = Mt) , (3.2) 

-d^ b (t) = -b 2 * 6 (t) , (3.3) 

and all functions satisfy boundary conditions (|3.ip . These functions form a complete system 
of (generalized) eigenfunctions: For functions v G L 2 (W + ) we have 

poo / roo \ 
V ^ = J ip b / p (tp)Tp b/p (sp)dp + y b (t)y b (s)\v(s)ds (3.4) 

in the sense of L 2 -convergence. This identity holds for continuous v G L 1 (M + ) n L 2 (R + ) 
and is extended first to L 1 (M + ) n L 2 (M + ) and then to L 2 (M + ) as in the case of the ordinary 
Fourier transform. 

We need the following technical result. 
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Lemma 3.4. For t G M + and b £ R we /iave 

Moreover, the function 

hit) = j\l -p>)«+W (^! C os(2tp) + -^sin(2^)) dp 

is uniformly bounded with respect to t > and b E M. satisfies 

r°° r°° I 1 if 6 = 

/ |I 6 (t)|(fc<C and / i|I 6 (t)|di<Cx { / x (3.5) 

Jo I ( 1 + ^ ) if o ^ 

u;zi/i C > depending only on the dimension. 

Proof. The first assertion follows directly from the definition of ipb since 

1 (l-6 2 )cos(2t) + 26sin(2t) 1 (1 - i&)V 2i + (1 + ib) 2 e~ i2t , , 

^ ) = 2 + 2(TT^ = 2 + i(TTO • (3 - 6) 

It is clear from the definition that lb is uniformly bounded. To establish decay in t we write 

M')=5/(i-p 2 r i,/2 fe^« i2 **. 

2 Jr P 2 + b 2 

and set G(p) = (1 - p 2 )^ +1)/2 and F 6 (p) = (p - ib) 2 /(p 2 + b 2 ). Let G and H b denote the 
inverse (distributional) Fourier transforms of G and Hb- 

It is well known that G(t) = CrfJ^/ 2+1 (|t|)|t| _a! / 2 ~ 1 , where Jd/2+1 denotes the Bessel func- 
tion of the first kind. The absolute value of this Bessel function behaves like t d l 2+l as t — > 0+ 
and is bounded by a constant times t -1 / 2 as t — > 00; see [AS64, (9.1.7) and (9.2.1)]. Hence, 
we have (^(t)! < Cmin{l, |i|~( d+3 )/ 2 }. Moreover, we compute that 
H b (t) = (27r) 1 / 2 5(t) - 2 3 / 2 ^ 1 /2| 6 | XR _ {bt) e -W . 

Thus we may rewrite Ib(t) in terms of G and Hb and get 



h(t) = \ ( G(2t - u)H b (u)du 
2 Jr 



1 /2 

G(2t) - (2vr) 1 / 2 |6| / G(2t - u)xr_(£w) e~^du 



1/2 f* 

|) G(2t) - (27T) 1 / 2 J G(2t + |) e^dn . 

In the last change of variables we have assumed that b ^ 0. From the bound |G(£)| 
\G{-t)\ < Cmin{l, |f|-( d + 3 )/ 2 } we easily derive that / °° \G(2t + u/6)|<ft < C. Moreover, 

1 r°° 

t\G(2t + u/b)\dt = - (t- u/b) \G{t)\dt 



4 



u/b 



s i(l |i||fi( ' )|4 + Bl |e( ' )|4 ) sc ( 1 + s)' 

This implies (|3.5p for b ^ 0. The case 6 = is similar. □ 



The next lemma establishes a connection between the function lb and the coefficient 
L { j \b) defined in IjTIIjl . 
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(2) 

Lemma 3.5. For L d '(b) we have the representations 

T m (h ,\C d J™h{t)dt for&>0, 

[Cd (Jo 00 h{t)dt + ir(b 2 + l)(^D/2) for b < . 1 ' ' 

(2) 

The function bt-tL),' (b) is countinuously differentiable. 



Proof. Because of the first bound in (13, 5p we may apply the dominated convergence theorem 
to write 

r I b {t)dt = lim [(1 - p 2 )( d + 1 )/ 2 r e ~ et2 ( P l - ^ cos(2tp)dt + -^-r sia(2tp)^ dtdp 
Jo 40 J J \p 2 + b 2 p 2 + b 2 J 



+ f(l-/) (d+1)/2 ^^(4)^P 



2p6 1 F /_p 

where F(x) = e~ x2 f Q x e y2 dy. Using the fact that 

lim — =F ' ^ 



e4-0 yfe \y/ej 2p 
we find 

for b / and / °° I 6 (i)cft = § for 6 = 0. By (fTTTT]) this yields d3T7j) . 

The fact that 6 i-> L ? (b) is C 1 away from b = is elementary. To prove continuity and 
differentiability at b = we again use dominated convergence together with the fact that 



lim /' 1 (l-P 2 ) (d+1)/2 7T^-7^ = ±- 



6->0±, 

We omit the details. □ 

3.3. Proof of Propositions [3TT1 and Lemma 13.21 After these preliminaries we can turn 
to the proof of local asymptotics on the half-space. We split the proof into three lemmas. 

Lemma 3.6. Under the conditions of Proposition EQ1 we have 

< 2C d [ I 4>\x){l - e d ) {d+m ^ 2 md (xdU/h) di d dxh~ d 
Jr^Jq 

+ irC d (b 2 + l)( d+1 )/ 2 f <f(x)*l /h (x d )dxh- d+1 - Tr (0tf+(6)<£)_ 

JR^ 

< Cl d ~ 2 h- d+2 (l + b d _T l min{6_, h/l}) , 
where C d is given in Here the constant C > depends only on d and Ca. 

Proof. First note that we may rescale <j) and thus assume I = 1 without changing the value 
of b. Since b is fixed throughout the proof we write H + instead of H + (b). 
To prove the lower bound we apply the variational principle and obtain 

-Tr(<£H+0)_ = Q mf <i Ti( 1 <t ) H + 4>) > q mf i (-Tr( 7 ^(F + )_ < /»)) = -T*(0(ff + )_0) . 
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Let a + {x,y) denote the integral kernel of (H + )_. From (|3.2p . (j3.3|) . and (I3.4|) we see that 
a+(x,y)=-i- J f m 2 -l)_e*'^-y'y h ^ d {xd d /h)^ 



and we get 

(rhH+rh\ < 

(2vr/i) 



Tr (</>#+</>)_ < — * / / </> 2 (*)(l£| 2 - 1)-^ (WO ^ 



+ " "*•+ 



1 

(2 7 r7i) d - 1 



Here we perform the ^'-integration and obtain the lower bound. 
We proceed to prove the upper bound. To simplify notation write 

f(x,0 = e ix '<'iP b/iidh) (x d C d ) 
F(x,H>) = e lx '<'* b/h (x d ). 

We define the operator 7 = (H + )°_ with kernel j(x,y) = 71 (z, y) +72(2,2/)) where 

4 

(2vr/i) rf y^ 6R d : |£|<i} 



71(2,2/) = 7— r-j / f(x,Z/h)f(v,Wdt, 



= (9 Md-i / F(x,?/h)F(y,Z'/h)d?. 

Thus, 7 satisfies < 7 < 1 and a variant of the variational principle, discussed in Appendix 
E yields 

- Tr(0F + <^)_ < Tr(</>7^+) = Tr(07i0F+) + Tr(07 2 <^H~ + ) . (3.8) 

We note that the range of (jyy4>, (frli^i an d 4'724 ) does not belong to the domain of H + . 
However, the functions 4>f and cj)F belong to the form domain i? 1 (M^_) of H + . Therefore 
(|3.8|) is valid if we interpret Tr((p'y(f)H + ) in the sense described in the appendix, namely 

Tr(</> 7l 0tf + ) = -i-j / qt\ct>f] # , (3.9) 

{2irh) a y {C6K d : | 5 |<i } 

where 

qtW\ = h 2 + hb\\4>{-M\l^-^l/i d (°) " H^II^( Ri ) , 

and similar for Tr(c/>72</>-ff + ). In the first summand we integrate by parts and use (|3.ip and 
(HU) to get 

\\n<t>f)\\l H ^ + ) = J Rd (f!^ 2 + IV0I 2 ) ^ 2 / fd (WO - \ ll^-.o)!^^ ^ d (0) . 

We insert this into (|3.9j) and due to (|2.4p and Lemma 13.41 we can estimate 
4 



Tr(^ 7l 0F+) < -— — j / / 2 (x)(|e| 2 -l)_^VtoM)^ + C/ i - d+2 . (3.10) 
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Note that the second summand in (|3.8f) is zero for b > 0. For b < we use (J3JJ) and (13.30 
to show that 



{2irh)( a L > J{|f|2<b2 +1} 



1 



/i 2 

+ 



— -— ^ / / |V^)|^ 7/l (x d )dx<. 

(27T/l)^ ^ J{|5'|2 <fe 2 + 1 } J R d_ I 

To estimate the last summand we use H^b/hH^ < —2bh~ 1 , H^b/ftHl = 1) an d (|2.4p to obtain 

/ |V^(x)| 2 ^ 2 /h (x d )dx < Cmin{-&//i, 1} . 

JR d + 

Performing the ^'-integration as before yields 

Tr(0 72 ^ + ) < - irC d h~ d+1 (b 2 + l)^' 2 f ^(x^l^dx 

JR^ 

+ Ch~ d+2 (l + bi^mmi^b-}) . (3.11) 



Here we also used the fact that 1 + (1 + b 2 )^ 1 ^ 2 min{6_, h} < C(l + b d S l min{6_, h}). 
Hence, the upper bound follows from (|3.8p . (|3.10p . and (|3.1ip . □ 

Lemma 3.7. Under the conditions of Proposition HOI we have 
2C d f f\ 2 {x){l-e d ) {d+l),2 ^ b/id {x d U/h)dUdx 

= L d 1] [ (t> 2 {x)dx + C d [ I b {t)dt [ ()) 2 (x',0)dx'h + r 1 {h,b) (3.12) 

JRf JO itf" 1 

with \n(h,b)\ < C(l + l/\b\)l d - 2 h 2 for 6/0 and |n(/i, 0) | < Cl d - 2 h 2 . For b < we a/so 

/ 2 (x)^ 2 /h (x (i )dx= / <p 2 (x',0)dx' + r 2 (h,b) (3.13) 

wrai/i \r2(h, b)\ < Cl d ~ 2 hbZ 1 ■ Here the constants C > depend only on d and C^. 
Proof. Recall that 

4 1] = T^w / (l^l 2 - = ^ Al - U) {d+1),2 dt d . 

(2vr) ct 7 R d Jo 

Hence, 

2C d [ f 4> 2 (x)(i - e d ) {d+1)/2 i> 2 b/id (xd d /h) d£ d dx 



R| Jo 

= 4 1} / </> 2 (^x + C d [ I </> 2 (x)(l - U) [d+l)/2 (2^1, (WO - l) di d dx . 
Jr* Jr^Jo v j 

We insert (j3.6|) and perform the £' integration and see that the right-hand side equals 

4 1} / ^(^)^ + C d [ <P 2 (x)I b dx , 
Jr± Jr± \ n / 



IN 
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with introduced in Lemma 13.41 To analyze the second term we insert 



d s cj) 2 {x' , s)ds 



(3.14) 



o 



and substitute Xd = th. We obtain 
2C d ' 



L 



(1) 



<a 2 (x)(i - e d ) {d+i)/2 ^ d (xd d /h) di d dx 

poo r 

<j) 2 {x)dx + C d h(t)dt 2 (x',O)dx'h 

Jo JRd- 1 



+ C 



oo rth 



JO 



d s (j) (x 1 , s)ds Ib(t)dtdx'h . 



Using (|2.4p and the remark at the end of Subsection 13. II we bound 

f [ d s <p 2 (x',s)dsdx' <Cl d - 2 ht. 

/R^- 1 Jo 

The first assertion of the lemma now follows from (|3.5p . 

The second assertion follows similarly by inserting ()3.14j) and by definition of D 

Note that the error terms in Lemma 13.71 diverge as b — > 0. Hence, we also need the 
following estimates that yield better results for |6| < Ch/l. 

Lemma 3.8. Under the conditions of Proposition [PI we have 



2C d 



<P 2 (x)(1 - ed) {d+l),2 ^l Kd (xdU/h) dUdx 



L 



(i) 



(j) 2 {x)dx + -L ( d l \ I (j) 2 (x',0)dx'h + ri(h,b) 



with \h(h,b)\ < Cl d ~ 2 h 2 (l + l 2 h- 2 \b\(l + | In |6||)). For b < we also have 

0< / <f) 2 (x)^ 2 b/h (x d )dx < C/ d /i~ 1 min{6_,/ir 1 }. 

Jr<*_ 

Here the constants C > depend only on d and Ga,. 

Proof. This proof is a variation of the previous one. Again, we write 



2C d 



wL Jo 



^(x)(l-^ 2 ) (d+1)/2 ^ d (^/^)^^=4 1) / cp 2 (x)dx 



+ C d [ <f> 2 (x)I b (^)dx. (3.15) 
Jw d v h J 



We add and subtract Jo to and from I b . According to the previous lemma and Lemma [3? 
we have 



Ca 



4> (x' , 0)dx'h 



< Cl d ~ 2 h 2 . 



Thus, it remains to control 
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Recalling the definitions of I b and Iq we see that the absolute value of this term is bounded 

by 

C [ <?{x)dx f\l-p^+^ b I±^ d p<Cl d \b\(l + \ln\b\\). 
jRd Jo p 2 +b 2 

This finishes the proof of the first assertion of the lemma. The second assertion follows 
similarly as at the end of the proof of Lemma 13,61 □ 

Proof of Proposition \3.1\ Combining Lemma [3761 with (13. 12ft . (|3.13p . and (|3.7p we obtain the 
first claim of Proposition 13.11 with a remainder 

\R hs {h,l,b)\ < Cl d - 2 h- d+2 (l + l&r 1 + (b 2 + i)(d+i)/2 b -i + b d r 1 mm{b-,hr 1 } 

< ci d - 2 h- d+2 \bi l {\ + \b\ + b d+1 ) . 

To obtain the second claim we combine Lemma 13.61 with Lemma 13.81 In this case the 
remainder is bounded by a constant times 

f-2 h -d+2 ^ + f h -2 m + | i n |6| |) + ((ft 2 + \)^)/H 2 h- 2 + b d _- 1 ) min{6_ , h/l]) . 

For |6| < h/l < 1 this simplifies to 

\R' hs {h,l,b)\ < Cl d - 2 h- d+2 (l + l 2 h~ 2 \b\(l + \ln\b\)) . 

This finishes the proof of the proposition. □ 

Proof of Lemma \3.2l Combining Lemma [3.61 with f|3. 15|) we obtain the claim with a remain- 
der bounded by 

■d+l 



\R'L(h,l,b)\ =C d <p 2 (x) iJ^f) dxh- d + <KC d (b 2 + iy d+i y 2 4?(x)* 2 b/h (x d )dxh- 
+ Cl d - 2 h- d+2 (l + b^ 1 min{6_, hr 1 } 



In the first term on the right side we substitute x d = th and use the first inequality in ()3.5 
to bound 

roo r 

cf) 2 (x',th)dx'\I b (t)\dt < Cl d - X . 







By Lemma 13.81 we also have 

< / (f> 2 (x)$> 2 b/h (x d )dx < Cl d K~ x min{6_, hr 1 } 

and the proof is complete. □ 

4. Local asymptotics close to the boundary 

Here we show how Proposition 12.31 and Lemma 12.41 follow from the results in Section [3j 
We straighten the boundary locally and estimate the operator H{b) given on Vt in terms of 
H + (b) given on the half-space MrL. 

In this section we work under the conditions of Proposition 12.31 Let <j) £ Co(^ rf ) be 
supported in a ball of radius I > and let inequalities (|1.10p and (|2.4p be satisfied. Then 
let B denote the open ball of radius I > 0, containing the support of 4>. Choose xq G B n dil 
and let u Xo be the inner normal unit vector at xq. We choose a Cartesian coordinate system 
such that xq = and u Xo = (0, . . . , 0, 1). 
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We now introduce new local coordinates near the boundary. Let D denote the projection 
of B on the hyperplane given by x d = 0. Since the boundary of fi is compact and in C , 
there is a constant Cq > 0, independent of xq £ d£l, such that for < / < Cq 1 we can find 
a real function / £ C 1 , given onDc R^ -1 , satisfying 

dntlB = {(x',x d ) : x £ D,x d = f(x')} f]B. 

The fact that d£l £ C 1 means that the functions V/ corresponding to different points xq 
and different values of I share a common modulus of continuity which we denote by u, that 
is, 

\Vf{x')-Vf{y')\<u{\x'-y'\) 

for all x',y' £ D. We assume that lo is non-decreasing and we emphasize that uj(5) | as 
510. 

The choice of coordinates implies /(0) = and V/(0) = 0. Hence, we can estimate 

sup \Vf{x')\ < supcj(|a/|) < w(Z). (4.1) 

x'eD x'&D 

We introduce new local coordinates given via a diffeomorphism p : DxM-> M. d . We 
set yj = (pj(x) = Xj for j = 1, . . . , d — 1 and y,i = </?d(x) = x d — f(x'). Note that the 
determinant of the Jacobian matrix of ip equals 1 and that the inverse of cp is defined on 
ramp = Dxl. In particular, we get 

p {dQ n B) C dM. d + = {y £ R d : y d = 0} . (4.2) 

Fix w £ ^(fi) with v = on M. d \ B. For y £ rant/3 put t^y) = v o yj~ 1 (y) and extend v 
by zero to R d . An explicit calculation shows that the effect of this change of coordinates on 
the gradient is small: 

Lemma 4.1. For v and v defined as above we have v £ H^(fi&i) and 



< Cuj{1) / \Vv(y)\ 2 dy. 



/ \Vv(x)\ 2 dx- / \Vv(y)\ 2 dy 
Jn Jm.^ 

Based on this estimate we now prove a result from which Proposition 12.31 follows. For 
£ Cg°(M d ) supported in B define <f> = </> o p 1 on ran p = D x R and extend it by zero to 
R d . It follows that 4> £ C^(R d ) and ||V^||oo < CI' 1 hold, with C depending only on C and 
uj. We set 6~ = \ni x ^dnnB 6(x) and 6 + = sup^g^^ b(x) and note that (6 + )_ < < o s , 

where b s was introduced in (j2.6f) . We also recall the notation H + (b^) introduced in Sectional 

Lemma 4.2. Under the conditions of Proposition [Ql i/iere is a constant > depending 
only on f2 suc/i i/iai /or < / < C^ 1 and < h < I we have 

Tr(4>H+(b+)4>)_ - Cl d h- d uj{l) (l + (6+) d+1 /ir 1 ) 

< Tr((f)H(b)<f))- 

< Tr{4>H + (b-)4>)^ + Cl d h~ d uj(l) (l + {b-f^hr 1 ^ . (4.3) 

Moreover, 

[ <j> 2 {x)dx = [ 4> 2 (y)dy, (4.4) 
Jn Jk^ 
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mi 



4> 2 (x)da(x) 



and 



no 



Lf{b{x))<j) 2 {x)da{x) -I%\b 



2 (y',O)dy' <Cl d - l u(lf, 
{2) ^ I 4> 2 (y',0)dy' 



(4.5) 



< Cl^ 1 ((1 + (6 ± ) d+1 )a;(/) 2 + (1 + (6 ± )l)/3(0) ■ 



(4.6) 



Proof. The definition of 4> and the fact that det J(p = 1 immediately give (|4.4p . In view of 
(|4.ip we can estimate 



(j) 2 (x)da(x) 



^ 2 (y',0)Vl + |V/|W< / 2 (j/,O)d l / + C'i« i - 1 a;(O 2 . 



This proves (|4~3|) . Using the fact that |L^ ) (6 ± )| < (7(1 + (6 ± ) d+1 ) we find 



4 2) (6(x))^(x)da(x)-L^(6 



P(y',o)dy' 



< 



Ml 



Lf(b{x)) - Lf(b±)\ <p 2 (x)da(x) + Cl d -Ml? (l + (& ± )- +1 ) 

The continuity of b, see (fTTOj) . and the fact that \-^L {2 \b)\ < C(l + bt) imply 

4 2) (6±) - Lf{b{x))\ < CP(l) (l + (6±)i) . 

Inserting this into the estimate above gives (|4.6p . 

To prove (|4.3p we first note that the variational principle implies 

Tr ((j)H(b + )(j))_ < Tr (<j>H(b)<j>)_ < Tr (4>H(b-)(/>) _ . 

Thus it remains to show that 

Tr tyHQ^W) _ -Tr(^ J ff + (6 ± )^)_ | < Cl d h- d u{l) (l + (6 ± ) d+1 /ir x ) 

To this end choose t> and u as in Lemma l4.1i First we estimate 

\v(x)\ 2 da(x)= [ |S(2/,0)|Vl + |V/| 2 <V> / \v(y',0)\ 2 dy' 



I an 

and using (|4.1|) 



f \v(x)\ 2 da{x) < {1 + Cu{lf) [ \v(y',0)\ 2 dy' . 
Jan jRd- 1 



(4.7) 
(4.8) 
(4.9) 



By decreasing, if necessary, the constant Cq from the beginning of this section we may now 
assume that I > is small enough such that 2Cu(V) < 1/2 holds. Then Lemma 14.11 ()4.4p . 
and (gi5]) imply, for 6 ± > 0, 



q b ±[v]>(l-Cuj(l))h 2 / |V?)(y)| 2 cfy + hb 

J9. d + 

= (l-2Coj(l))q+ ± [v) 

+ 2Cu(l) ( ^ / |Vv(y)| 2 dy + /i6 ± / 
= (l-2C W (Z))g+ ± [5] + 2Ca;(Z)g+[T;], 



\v(y',0)\ 2 dy' 



v(y',0)\ 2 dy' 



\v(y)\ 2 dy 



\v(y)\ 2 dy 



(4.10) 
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where q + is the same form as q + but with h replaced by h/y/2. For < we get, using 



Q 6± H>(1-Ca;(0)^ / |V5(y)|^y 

+ {l + Cu(lf)hb ± I \v(y',0)\ 2 dy'- I \v(y)\ 2 dy 



>(1 - 2Cw(/))g+ [S] + 2CMZ)g+ b± [5] . (4.11) 
To deduce estimates for Tr (<f)H(b )(/>)_ we recall the variational principle 
-Tr UH{b ± )(j))_ = inf Tr (070F(6 ± )) , 

where we can assume that the infimum is taken over trial density matrices 7 supported in 
B x B. Fix such a 7. For y and 2 from DxR set 

so that < 7 < 1 holds. Moreover, the range of 7 belongs to the form domain of cf>H + (b^)4>. 
First, we assume 6^ < 0. According to (|4.1ip it follows that 

Tr ((/r/^HQ^)) > Tr (faf ((1 - 2Cu{l))H + {b±) + 2Clo{1)H + (Cb^) 

> - (1 - 2Cw(/))Tr f0i/+(6 ± )0 N ) - 2Cu(1)Tt (ipH+iCb*) 



where the operator i? + is generated by the form q + . This implies 

Tr(0^(6 ± )0)_ < Tr(^+(6 ± )0)_ + 2Cw(0Tr (^H+iCb*) 
and Corollary 13.31 yields 

Tr(^(6 ± )0)_ < Tr(^i/ + (6 ± )0)_ + Cl d h- d Lo(l) (l + (b±) d+1 h/l 
for 6± < 0. 

In the same way we can treat non-negative b^ using (|4.10p and we obtain the lower bound 
in (14.71) . Finally, by interchanging the roles of H{b^) and H + (b^), we get an analogous 
upper bound and the proof of Lemma 14.21 is complete. □ 



Proof of Proposition \2.3\ and Lemma \2.J\ The assertions follow from Lemma 14.21 together 



with Proposition 13.11 □ 

If we combine the estimates of Proposition 12.21 Corollary 13.31 and Lemma 14.21 we obtain 
the following simple bound that is useful to estimate error terms. 

Corollary 4.3. There is a constant Cn > with the following property. Let (ft £ Cq° be 
supported in a ball of radius I > and let (12.41) be satisfied. Assume that b is a real constant 
independent of x . 

Then for < I < Cq and < h < I the estimate 



Tr (<t>H(b)<j>)_ < Cl d h~ d (l + b^hV 1 
holds with a constant C > depending only on d, and u. 
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5. Localization 

In this section we construct the family of localization functions (4>u)u&R d an d prove Propo- 
sition 12.11 The key idea is to choose the localization depending on the distance to the 
complement of Q, see [H6r85, Theorem 17.1.3] and [SS03 for a continuous version of this 
method. 

Fix a real-valued function cf) G C^°(M. d ) with support in {|x| < 1} and \\4>\\2 = 1- For 
a,x £ R 11 let J(x, u) be the Jacobian of the map u \— > (x — u)/l(u). We define 



such that (f) u is supported in {x : \x — u\ < l(u)}. By definition, the function l(u) is smooth 
and satisfies < l(u) < 1/2 and HVZH^ < 1/2. Therefore, according to [SS03] . the functions 
4> u satisfy flOJ) and <^ for all u G M d . 

To prove the upper bound in Proposition 12. 1|, put 

7 = / ^ (<^iT(6)<^)° <l> u l{u)- d du. 



Obviously, 7 > holds and in view of (|2.3p also 7 < 1, hence, by a variant of the variational 
principle discussed in the appendix, 

-Tr(H(6))_ < Tr (7/7(6)) = - / Tr (<j> u H(b)<j> u )_ l{u)- d du . 

To prove the lower bound we use the IMS-formula. For <fi G Cg°(M d ) and v G i7 1 (Jl) we 
have 

^Vv • V (0 2 U) + ^Vv ■ V (<j) 2 v) = |V(0v)| 2 - \V<f>\ 2 \v\ 2 . 
Combining this identity with the partition of unity (|2.3|) yields 

(%[<M - {v,h 2 (V<j) u ) 2 v) L2{n ^ l{u)- d du. (5.1) 

Using (|2.2p and (j2.3|) one can show [SS03], for every x G 

(V0 u ) 2 (x)^(u)^(iu < C f ^ 2 u {x)l{u)~ d - 2 du. 

We insert this into (|5.ip and deduce 

Tr(H(b))_< [ Tr(<j) u (H(b)-Ch 2 l(u)- 2 )<f) u )_l(u)- d du, (5.2) 
Jn* 

where Q* = {u G M. d : supp</> u HQ 7^ 0}. For any u G R, let /0 U be another parameter 
< p u < 1 and estimate 

Tr {MH(b) - Ch 2 l(u)- 2 ) ( f )u )_ < Tr + Tr (</>„(/>„# (6) - Ch 2 l{u)- 2 )<t> u ) _ . 

We now claim that choosing p M proportional to h 2 l(u)~ 2 yields 

Tr (^(g(b) - Ch 2 l(u)- 2 )<f> u )_ < Tr + g-pz^ ( 1 + 7 ~) J • ( 5 - 3 ) 

To see this, let us write t u = p u /(p u + Ch 2 l(u)~ 2 ) and note that r u < 1 and 

Tr {MpuH(b) - Ch 2 l(u)- 2 )<p u )_ = Ch 2 l(u)- 2 (l - r u )- 1 Tr(^F(V^6)</) u )- . 
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Here H is generated by the same quadratic form as H but with h replaced by s/r^h. If 
4> u n dVt / 0, we have Iq/A < l(u) < lo/V~3, see (|2.1U|) and (|2.11|) . and we can apply 
Corollary 14.31 to estimate 

Tr(<p u H(^b)<p u )- < Cl(u) d h- d T- d / 2 (l + {b m ) d + l hl{u)- 1 ) . 

With our choice of p u proportional to h?l(u)~ 2 we find that r u is order one and (|5.3p follows. 
If (p u E Cq° we can argue similarly by using the lower bound in Proposition 12.21 and get 

Tr {K(H(b) - Ch 2 l(ur 2 )ct>u) _ < Tr (cj )u H{b)ct> u )_ + . (5.4) 

Finally, we insert ()5.3[) and (|5.4p into (|5.2p and arrive at 

Tr (H(b))_ < [ TT{(t) u H{b)(j) u )_l{u)- d du + Ch- d+2 [ l( U y 2 du 
Jn* Jn\u 

+ Ch- d+2 J (l(u)- 2 + (b m ) d _ +1 hl(u)- 3 ^ du , 

where U = {u£R d : dQ.C\B u ^ 0}. Thus the claim of Proposition O follows from (pT[5|) 
and (I2TT6]) . 

Appendix A. A geometric lemma 

In the proofs of Theorem 11.11 and Theorem 11.31 we used the following estimate. 

Lemma A.l. For every domain O C M. d with dQ E C l there is a constant C with the 
following property. For every < Iq < 1 and u E M. d let l(u) be defined as in (|2.ip fry 

1 / — 1/2\ — ^ 

|(«) = - M + (dist(«,R d \fi) 2 + Z^ J . 

T/ien /or any relatively open N C (90 i/ie set 

U* = \u E R d : dist(u, 90) < l(u) A dist(u, 50 \ N) > Z(u)} 

safc/ies 

limsup-|?7*| d < Co-(iV). 

Z 4-0 '0 

ifere | • \d denotes the d- dimensional Lebesgue measure on lR d and o~(-) denotes the d — 1- 
dimensional surface measure on dVt. 

Proof. We split £T into two parts [/* = U* n and ?7* = f7* n M d \ and we prove the 
assertion separately for each of them. We begin with U* . Note that for u E we have 
dist(ii,M rf \ 0) = dist(u, 50). We first argue that there is a constant Li such that 

[/* = {u E : dist(u, dU) < L lo A dist(u, dU\N)> l(u)} (A.l) 

and such that Iq/4 < Li < Zo/v^- 

To prove (jA.ip let us consider the function 

*io (*) = J f 1 + + *o) _1/2 ) _1 " ^ > x > . 
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This function is continuously differentiable and satisfies Fi (Q) = /o/(2(Zo+l)) > 0, Fi (x) < 
for x > 1/2, and 

K (*) = f (* 2 + «) ~ 1/2 (i + (* 2 + ^o 2 ) 1/2 ) " 2 - 1 < -\ 

for all x > 0. Hence, there is a unique L; G (0,1/2] with Fi (Li ) = 0. Moreover, since 
F io (/ /4) < < F, (Zo/V3), we have Z /4 < L lo < l /V3. 

By definition, all u G O with dist(«, 90) = L; satisfy F/ (dist(n, 00)) = 0, thus l(u) = 
dist(/u, 00) = Li . The fact that Fi is decreasing shows that the inequality dist(u, 00) < L[ 
implies i 7 ) (dist(u, Oft)) > 0, thus dist(u, 0O) < l(u). Similarly, the inequality dist(«, 00) < 
l(u) implies dist(it, 00) < Li . This proves (jA.lj) . 

Our next step is to fix an < e < 1 and to decompose U* = Ut, U U* with 

C/> = {u G O : dist(«,0O) < (1 - e)L lQ A dist(u,0O \ JV) > 

U* = {u G O : (1 - e)L < dist(u, 50) < L, A dist(n, 00 \ iV) > l(u)} . 

Thus, 

\U*\ d <\U^\ d + \U:\ d . 
The second term on the right side can easily be bounded, 

\U*\ d < \{u G O : (1 -e)L < dist(u,0O) < L}\ d < f '° <r (0O*) eft < C/ e • 

Here we wrote 0O( = {a £ SI : dist(u, 0O) = i] and used the facts that a{dVtt) is uniformly 
bounded and that L\ Q < Iq/ \/3- 

After these steps we have reduced the lemma to proving that 

limsup -\Ul\ d <Ca(N) (A.2) 
i io h 

with a constant C independent of e. To do so we start from the representation 

Wt\d = / <r(Ut)dt, (A.3) 

Jo 

where 

[/;* = {iiell : dist(u, 00) = t A dist(tt, 0O \ N) > l(u)} , < t < (1 - e)L io . 

Recall that every u £ U* and, in particular, every u G £/> satisfies dist(it, 0O) < Z(u). We 
now claim that for every < e < 1 and every < Iq < 1 there is an r > such that every 
u G £/> satisfies 

l(u) > dist(u,0O) + r. 

This follows again from the monotonicity and continuity of the function Fi . Indeed, we can 
setr = F /o ((l-e)L /o ). 
We consider the set 

N := \J \J \J {y G R d : \y - x\ < r} n 00 

0<t<(l-e)L; nSC/ t * xedn, |a;-tt|=t 

and show that 

N C N (A.4) 
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and 

a(dN) = 0. (A.5) 

To prove (|A,4h let < t < (1 — e)L/ , x,y £ dfl with |x — y| < r and u G J7 t * with 
\x — u\ = t. Then 

\y — u\ <\y — x\ + \x — u\ < r + dist(ii, dVL) < l(u) . 

Since dist(u, d£l \ N) > l(u) by the definition of U*, we infer that y G N. This proves (|A,4p , 
To prove (|A.5j) we note that N satisfies the following uniform interior ball condition. For 
each y G dN there is an open ball B C M. d of radius r such that y G dB and -B n d£l C iV. In 
order to prove (|A.5P we introduce local coordinates similarly as in Section HJ In this way we 
are reduced to the situation where N is a subset of satisfying a uniform interior ball 
condition (with a possibly smaller radius). The claim (|A.5P follows from Lemma I A . 2 1 b elow . 
The definition of N easily implies that 

[/* c £/" t * := G O : dist(«,<9fJ) = f A dist(n,iV) = t| 
for all < t < (1 — e)Li Q . Moreover, we can estimate with a constant depending only on 
a(Ut) < C (a(N) + a({x G dQ, \ N : dist(x,A>) < t})) 

< C (a(N) + a({x G d£l \ N : dist(x, N) < / })) . 
The second bound used (|A.4h as well as (1 — e)Li < (1 — e)/o/v / 3 < Zo- Thus, from ()A.3[) . 

|f/>U < Cl (a(N) + a({x G 50 \ N : dist(x,A>) < l }j) . 
Therefore, in order to prove ()A.2p . it remains to estimate 

a({x eOQ\N: dist(x, N) < l }) = f X i Q (v)dv(v) > 

Jan 

where x; denotes the characteristic function of {x G 30 \ N : dist(x, N) < Iq}. We note 
that lim; ^o Xi = Xqn pointwise. Thus, the dominated convergence theorem and (1A.5|) imply 
that 

lim a({x G dO, \ N : dist(x, N) < l }) = . 
loiO 

This completes the proof of (|A.2h , 

For U* we get an analoguous bound by following the same strategy. In this case the 
estimates are somewhat simpler since, for u G M rf \ O, we have l{u) = + 1) and this 

plays the role of Li . □ 

Lemma A. 2. Let A C W 1 be bounded. Assume that there is p > such that for each x G dA 
there is a ball B C W l of radius p with x G dB and B C A. Then \dA\ n = 0. 

Proof. Let 5 > be a constant to be specified later and put l m = 5p5~ m for m > 0. We 
denote by Q m the collection of open cubes of side length l m centered at points in (l m Z) n . 
Let C m be the collection of those cubes in Q m that intersect both A and M n \ A. Since A 
is bounded, v m := #C m is finite. We claim that for all sufficiently small 5 > there is a 
constant M < 5 n such that for all m > 1 



(A.6) 
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Deferring the proof of this bound for the moment we now explain why it implies the lemma. 
First, we iterate (|A.6[) to learn that v m < M m UQ. Thus, since dA C UoeC m Q ^ or an y m ' we 
conclude that 

\dA\ n < \Q\n = l>m < 5 n p n (5- n M) m is ^ as m y oo . 
QeCm 

This proves |cL4| n = and we are left with showing (|A.6[) . 

To do so, we fix m > 1 and an arbitrary cube Q C C m —\. When passing from m — 1 to 
m, this cube is subdivided into 5 n cubes in Q m . We shall show that if 5 > is sufficently 
small then at least one of these cubes of side length l m does not belong to C m (i.e., does not 
intersect both A and W 1 \ A). This will imply with M = 5" - 1. 

Consider the cube Q' G Q m in the center of Q. If this cube does not belong to C m we 
are done. Thus, we may assume that Q' intersects both A and W 1 \ A. Because of our 
assumption on dA there is an open ball B of radius p such that B C A and dB HQ' ^ 0. 
We now make use of the following 

Claim. There is a constant C n > such that if B C M. n is an open ball of radius r > C n 
with BnQ / 0, where Q = (-1/2, l/2) n , then 7 + Q C S for some 7 G Z n with < 2. 

Indeed, one can take C n = maxj-^/n, n/2}. The proof of this claim uses only elementary 
geometric facts and is omitted. 

By a rescaled version of the claim we infer that, under the assumption that p > C n l m , 
there is a cube which is contained in B and whose center is at most an oo-distance 2l m away 
from that of Q'. Since Q' lies in the center of Q this cube is also contained in Q. Moreover, 
since it is contained in B, it is also contained in A and, therefore, does not belong to C m . 

Finally, we argue that for all 5 > small enough the assumption p > C n l m is satisfied for 
all m > 1. Indeed, this assumption is equivalent to 1 > C n 55~ m , which holds uniformly in 
m > 1 provided we choose 5 < 5C" 1 . This completes the proof. □ 



Appendix B. A variant of the variational principle and a sharp bound on 

Tr(-A 6 - A)_ 

Here we mention the following extension of the variational principle that we used in the 
proof of Proposition 13.11 

Let (M,p) be a measure space and let (f a )aeM be a measurable family of functions in a 
separable Hilbert space Q, such that 

\(iPJ a )\ 2 dp(a)<\\iP\\ 2 (B.l) 

M 

for all ip G Q. Assume that A is a self-adjoint, lower semibounded operator in Q with 
quadratic form a such that 

f a £ dom[a] (B.2) 

for all a £ M. 

Let the operator 7 in Q be given by 7^ = J M (f a ,ip)fadpi(a). Then 7 satisfies < 7 < 1. 
Let us introduce the notation 

Tr^7 = a [f a ] dp(a) . 
Jm 
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Then we have 

- Trvl_ < TrA-y , (B.3) 

provided J M a[f a ]-dfi(a) < oo. 

Let us illustrate these notions by adding the following sharp estimate, a simple form of 
the upper in Proposition 13.11 which is based on a method introduced in [Kr692] . Here we 
only assume that the boundary of C M. d is Lipschitz continuous and that — A& is generated 
by the quadratic form given in (jl.ip . 



Proposition B.l. For <f> G C^(R d ) and A > 

Tr(0(-A 6 -A)0)_ >4 1} A 1+d / 2 f \(j)(x)\ 2 dx 



(2vr) c 



A d / 2 ( / b{x)\<P(x)\ 2 da{x)+ / \V0\ 2 dx 



on 



Proof. To adopt the notation introduced above, we set Q = L 2 (Q), M = {£ G R d : |£| 2 < A} 
and ii to be Lebesgue measure. If we choose f^(x) = (2ir)~ d ^ 2 e lx '^ then (jB.ip and (|B.2|) are 
satisfied and the claim follows from (jB.3|) . □ 

If we choose <j) = 1 on O, we get 

Tr(-A 6 -A)_>^|n|A 1 +-/ a - 7 |^ / 6(x)Mx)A d / 2 . 

Jan 

This generalizes the bound proved in [Kro92] for the case of Neumann boundary conditions. 
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